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Abstract 

We study the geometry of a semiflexible polymer at fi- 
nite temperatures. The writhe can be calculated from 
the properties of Gaussian random walks on the sphere. 
We calculate static and dynamic writhe correlation func- 
tions. The writhe of a polymer is analogous to geometric 
or Berry phases studied in optics and wave mechanics. 
Our results can be applied to confocal microscopy studies 
of stiff filaments and to simulations of short DNA loops. 

1 Introduction 

Experiments which micro-manipulate stiff polymers, 
such as actin filaments or DNA have lead to a re- 
newed interest in the static and dynamic proper- 
ties of semiflexible polymers. This article treats the 
writhing properties of stiff polymers, in particular 
those which have a length comparable to the persis- 
tence length. The static and dynamic writhing prop- 
erties discussed here can be easily measured using 
confocal microscopy on fluorescently marked actin fil- 
aments for which the persistence length is 15|j.ra, and 
which are readily available in a range of length vary- 
ing between 5|j.ra to SO^ita. 

In this article I concentrate almost exclusively on 
the writhing fluctuations of polymers; clearly a stiff 
polymer has both bending and twisting degrees of 
freedom, but it is known y, || || that the relaxation 
times of the twist modes are much faster than those 
of bending modes, which drive the writhe dynamics. 
Thus the longtime effective degrees of freedom of a 
fluctuating stiff polymer are those of writhe discussed 
here. 

At finite temperatures the trajectory of a semiflex- 
ible filament is not smooth and direct application of 



the standard relationships between, for instance geo- 
metric torsion and writhe turns out to be quite sub- 
tle. Indeed the usual definition of torsion requires 
curves that are C^, whereas the trajectory of a ther- 
malized stiff polymer is much less smooth than this; 
in fact both the torsion and its integral are divergent 
so that the usual relationship between writhe and tor- 
sion, InWr + J Tds — mod In, is useless. Here we 
explore representations of the trajectory of a polymer 
based on random walks on the sphere which have the 
advantage of being easily visualized. We are able to 
use this representation to calculate the static and dy- 
namic fluctuations of a stiff polymer. We discuss the 
analogy between Fuller's result for the writhe of a 
polymer and the geometric phase known from quan- 
tum and optical physics. 

The main quantitative predictions made in this 
paper are for the static and dynamic correlation 
functions for the writhe, Wr in several different 
geometries. In particular we calculate (Wr^) and 
((yyr(t) — yVr(O))^) for open and closed semiflexible 
polymers. With these results we are able to under- 
stand quantitatively a number of results for writhe 
autocorrelation functions which have been measured 
recently in numerical models of DNA Q. We also 
show that short, highly curved sections of polymer 
give enhanced contributions to the writhe of a poly- 
mer. We explicitly calculate the probability distri- 
bution for writhing, Vy, ( Wr) for short polymers and 
show that it is strongly non-Gaussian. 



2 Writhe of a curve 

Consider fig. ™). It represents a solid bar that has 
been bent three times by 90°. On one face of the 
bar I have added an unit arrow, n, which is perpen- 



3 GEOMETRY OF THREE DIMENSIONAL CURVES 



dicular to the local tangent of the bar, t. Following 
the motion of n we notice that after the triple bend 
there has been a rotation of the arrow by 90° about 
the vertical axis despite the fact that t has come back 
to its original direction. This rotation is due to the 
writhe of the path. Mathematicians measure writhe 
in units of 2n. Thus fig. (|l|) corresponds to a writhe, 
Wr of 1/4. 



dimensional surface this angle is given by 




Figure 1: A bar is bent three times, leading to a rotation 
of 90° of the vector n marking one face of the prism. 
This rotation is due to the writhing of the shape in three 
dimensions. 

The writhing in fig. (|^) can be best understood 
by studying the geometry of the tangent vector t to- 
gether with the normal vector n. Since the vector t 
has unit magnitude, is thus lives on a unit sphere, Sz, 
drawn in figure (||) . The vector n is perpendicular to 
t: it must always remain tangential to the sphere. 
The geometry of the bent bar in fig. (p has been 
translated onto the sphere in fig. (||). Start at bot- 
tom right of fig. (|j) and move vertically, following the 
vector t. The whole vertical section, A corresponds 
to a single point on the sphere, the north pole. Pass- 
ing through the first bend a corresponds to the arc, 
also labeled a on the sphere. The second straight sec- 
tion of the bar, B points to the left and corresponds 
to the point B on the left most edge of the sphere. 
Despite the discontinuities in the direction of bar the 
vector n is always continuous on Sz ■ The vector n is 
parallel transported around a loop on the sphere. 

The rotation of the vector n is a manifestation of 
the curvature of the sphere. Parallel transport of a 
vector around a closed path on a surface leads to a 
angle between the initial and final vectors. For a two 



Aa = 



JC dS raodulo 2n , 



(1) 



where /C is the Gaussian curvature of the surface and 
the integral is over the region enclosed by the path. 
For a sphere the Gaussian curvature is equal to unity 
so that the angle of rotation between two ends of a 
bent prism is calculated from the area enclosed by 
the trajectory on the sphere. This is just the con- 
tents of Fuller's result [m. In our simple example the 
solid angle enclosed by the the curve (a, P, y) is n/2 
leading to the rotation of n by 90° . 

Any real prism, made of an elastic material has 
an additional degree of freedom associated with it: 
At any point it can be twisted, at angular velocity 
cutls) about the tangent, compared with the parallel 
transported frame shown in figures (|l|,||). Clearly the 
total rotation in the laboratory frame is given as the 
sum of the rotation due to the writhing, plus the 
rotation J cut ds coming from the relative motion of 
the parallel transported and twisting frame. This is 
the content of White's theorem Q for a curve. Cn — 
Wr -|- Tw, where 2nCn is the total rotation, IttTw 
the rotation due to twisting and 2n'Wr the rotation 
due to writhing. 



3 Geometry of three dimen- 
sional curves 

In this section I shall review some of the elementary 
properties of a line in three dimensions, pointing out 
the relationships between different representations of 
the geometry. In particular I will discuss the classic 
choice of the Frenet frame and compare it with the 
parallel transported frame of fig. (1,2) to show how 
and when the two descriptions can be related. 

At each point of the curve associate three orthog- 
onal unit vectors {t,n,b,}. The vector t is tangent 
to the curve, while n and b are perpendicular to the 
tangent, they both live in the surface of the sphere of 
fig. (0) . A theorem of Euler H then shows that the 
motion of this frame as a function of s, the curvilin- 
ear distance, is a pure rotation. Thus grouping the 
vectors into a matrix, M gives 



dM 
~ds" 



Q.M 



(2) 



AD 




Figure 2: Trajectory in the tangent space. The points A, 
B, C, D correspond to the straight sections of the prism 
in fig. (hi). The three bends correspond to the geodesic 
arcs a, |3, y. The vector n is parallel transported along 
the path, leading to an angle of 7t/2 between the initial 
and final vectors n. The vector n' (corresponding to 
the Frenet frame) always points in the direction of 
motion on the sphere. 

with D. anti-symmetric. Different choices are possi- 
ble; we shall now discuss two: firstly the classic Frenet 
frame, known from elementary treatments of line ge- 
ometry in three dimensions, secondly the choice of 
parallel transport corresponding to fig. ([l|). 

The most common choice for the rotation matrix, 
is the following corresponding to the Frenet frame m , 



(3) 



This representation is particularly useful because 1/k 
has a simple geometric interpretation: it is the radius 
of curvature of the line at the point s. The element 
T is known as the torsion. 

As shown above the geometry of three dimensional 
curves becomes easier to understand if we work in the 
tangential space Sz- What is the geometric picture 
that corresponds to the choice of the matrix in eq. 
(^)? On the sphere, the curve t(s) moves at speed k 
in the direction n' which is parallel to dt/ds. As a 




function of s the vector n' rotates in the surface at 
angular velocity (x>t = t. Since we are particularly 
interested in the internal geometry of the sphere let 
us now transform the parameterization of the curve 
from unit speed in real space to to a unit speed curve 
on the sphere. The element of length transforms as 
dcT — Kds; the angular velocity of n' expressed as a 
function of cr (rather than s) is t/k, j^. In fig. (||) 
t/k is zero along the arcs a, P, y but at the corners, 
A, B and C there is a rapid rotation of n' by 90°, 
corresponding to a delta-function singularity in the 
angular velocity. The representation becomes singu- 
lar in for curves which are not sufficiently smooth, or 
for curves where k passes through zero. 

The second representation of the geometry (some- 
times known as the Fermi- Walker frame), |l^ is the 
following 

/ cc |3 \ 
n = -a . (4) 

V -P / 
It corresponds to the geometry shown in figures (0, 
H). In this frame the vectors n and b are parallel 
transported: The original definition |ll| of parallel 
transport of a vector v (due to Levi-Civita) is that 
the rate of change of the vector, projected back on to 
the tangent plane is zero. For our space curve this 
corresponds to the equation 



dv 
da 



t t.- 



dv 
'dCT 



0. 



(5) 



We see that both n and b (as well as an arbitrary, 
constant linear combination of b and n) evolving with 
the matrix eq. (Q) obey this equation. In figure (g) 
the vector n is translated in such a way that dn/dcr 
is always normal to the surface, in agreement with 
this definition. 

For sufficiently smooth curves both representations 
are equivalent. At each point on the curve perform a 
rotation about the vector t in the following manner: 



n' 
b' 



cosO 
— sin( 



sm t 
cost 



(6) 



where n' and b' are the vectors in the description 
eq. (^ and n and b are the vector in the description 
of eq. (||). The two representations are be linked by 
setting 

e = T (7) 

a — KCOS0 (8) 

P = Ksine . (9) 



4 EXAMPLES OF WRITHE 



Consider an arbitrary trajectory on the sphere 
comparing the two frames, starting from a common 
initial vector n(0). After the loop there has been a 
rotation of 



T ds 



dCT 



(10) 



between the two description. If the angle between 
the parallel transported frames n(0) and n(L) is Aq 
then the angle between n(0) and n'(L) is (Ao + 9l)- 
Clearly this must be equal to a multiple of Ztt since 
n'(0) and n'{L) are identical (since they are both 
equal to dt/dff). Thus we deduce 



Aq + 



dcr = 2nn , 



(11) 



which when combined with equation (Q) is the Gauss- 
Bonnet theorem for a sphere. 



4 Examples of Writhe 

4.1 Bent Bar 

The trajectory on the sphere, shown in in fig. (0) has 
three rotations of n/2 corresponding to each bend in 
the bar thus J{t/k) dcr — 3n/2. Aq is also equal to 
n/2 as is the area enclosed by the curve on the sphere 
as required by eq (|l|) . Thus 2nWr + J t ds = 2n. 

4.2 Circular Helix 

A second example is a long helix, parameterized by 
the equation 



r u 



(acosu, asmu, uj 



for the case a small. The torsion is given by 

1 



1 +q2 



(12) 



(13) 



with ds — duV 1 + Q'^ . As a function of s the trajec- 
tory of t on the sphere is a circle or radius a. Con- 
sider the various contributions to eq. (11 ) for a line of 
length L — IttN vT+c? , so that the are N complete 
helical repeats. The integral of the torsion is given 

by 



ds 



2N7T-N7ta^-FO(a'*l 



(14) 



The result is at first very surprising. In the limit 
a — > the helix becomes a straight line, however the 
torsion and its integral remains finite: The straight 
reference state is singular for the calculation of the 
torsion of a curve. 

This second example shows that the torsion is sen- 
sitive to small scale details in the representation of a 
curve, in particular to small scale helical structure (as 
is found in stiff biomolecules) . At finite temperature 
the situation is even more delicate. Even for a fila- 
ment with a linear ground state small scale structure 
is excited by thermal fluctuations so that the integral 
of the torsion is badly behaved. In fact the stan- 
dard relationships between writhe and torsion, lead 
to large, uninteresting zero order contributions which 
mask the interesting higher order term, as found here 
for a simple helix. 

The enclosed area by the path t(s), is better be- 
haved than the torsion so that we can use eq. (P 
and the fact that the circle of area ttq^ is traversed N 
times to find Aq = Nttu^. The total rotation, which 
is physically interesting, corresponds to the second 
order term in a^ appearing in eq. (|l^) . We conclude 
that it is better to calculate the writhe from Fuller's 
result on the enclosed area of the curve t(s), rather 
than from the integrated torsion. We now generalize 
these remarks to finite temperature. 

4.3 Local Geometry 

Locally a smooth writhing curve can be written in 
the following form 



r u 



KT , K -, 

-^s,-s,s 



(15) 



showing that kt measures the rate at which the curve 
becomes non-planar. We have chosen local coordi- 
nates so that the tangent points in the S^ direction. 
From this expression we again see that a straight fila- 
ment with K = can be approached with an arbitrary 
value of the torsion t. 

On ^2 we find that t « ( — kts^/2, ks, 1 ). When we 
transform to unit speed on 52 using ks = cr we find 



^ 2 1 
- — o-^,o-, 1 

2k 



(16) 



so that the projection of the path from the sphere 
onto the tangent plane is curved, with radius of cur- 
vature k/t. Thus t/k has a natural interpretation on 
the sphere and is indeed the geodesic curvature of the 
trajectory on 52. 



5.2 Short Polymers 



5 Bending of stiff polymers 
5.1 Energetics 



The bending energy of a stiff polymer is written in 
the form 



-f 



ds 



ds 



where K is the bending modulus of the polymer. Let 
the energy be measured in units of the temperature 
then K is the persistence length of the polymer. This 
energy is completely analogous to that used in the 
theory of flexible polymers M% 



-Hex — -T 



/dr(s] 
V ds 



ds 



(18) 



except that it is written in terms of the tangent vector 
t rather than the real space position r. The energy 
( p!8| ) describes a Gaussian polymer with end to end 
distance (R^) — 3L/L. Like a flexible polymer the 
distribution of density P(t,s) obeys a diffusion like 
Fokker-Planck equation [O . 



9P(t,s) 1 7 , . 



(19) 



where Vj. is the Laplacian operator on the sphere. 
As a function of s the vector t diffuses with diffusion 
coefficient 1/(2K). The typical separation on 52 of 
two points separated by a curvilinear distance s in 
real space is 

r2 = 2s/K, (20) 

when s <C K. 

From the result eq. (19) we understand why tor- 
sion, T is poorly defined at finite temperatures: Con- 
sider a discretization of a semiflexible polymer with 
sections of length a in real space. On the sphere the 
polymer is represented as a random walk with points 
separated by a distance comparable to \/a/K. The 
direction of each step on the sphere is completely un- 
correlated with its neighbors. Thus the angle between 
successive links is of order of one radian. There are 
L/a sections so that the torsion locally has a value 
T ~ ±l/a and jTds ~ ii/L/a which diverges with 
the cutoff. The physically interesting writhe is thus 
hidden by the diverging zero order fluctuations. This 
is similar to the example of the helix above where the 
physically interesting writhe was hidden by the small 
amplitude helical pitch. We shall now show that de- 
spite some subtleties the angle between the vectors 



n(0) and n(L) can be calculated from the area en- 
closed by the random walk on Sz- 

This result can be linked with the natural class 
of functions for which the Frenet frame is defined. 
Bishop ||l^ showed that the Frenet frame (and thus 
the torsion) requires a function space of C curves. 
The frame defined by parallel transport has lower 
smoothness requirements. Clearly a semiflexible 
polymer, which is only C^' , is too rough for the 
Frenet frame to be useful. 

5.2 Short Polymers 

Consider a polymer with the following boundary con- 
ditions: The ends of the filament r(0) and r(L) are 
free to move but the chain is constrained by external 
couple so that t(0) and t(L) point in the ^ direction, 
as in fig. (|l|). The trajectory is closed on the sphere 
(but open in real space). If the polymer is short, 
(L <ti K) the trajectory never moves too far from the 
north pole; locally the geometry of the sphere remains 
Euclidean. The writhe of the filament can thus be 
calculated from the area enclosed by a closed Gaus- 
sian loop in a plane. This problem has been treated 
by Levy Q. The enclosed area, fig. (Q), scales in 
the same manner as Rj. but is randomly signed (de- 
pending on whether the path is traversed clockwise 
or anti-clockwise). The mean squared area has a non- 
zero average: 



(Wr^ 



48K2 



(21) 



In fact the whole probability distribution is known 
and is given by 



PkOVt) 



7tK 



1 



L cosh^(27rWrK/L) 



(22) 



Although the distribution of the size of the loops is 
Gaussian the distribution of the writhe has wings 
which are simple exponentials, implying that the free 
energy of a strongly writhing state varies linearly with 
Wr, 

F = 47tK Wr/L . (23) 

This result implies that there is a critical torque 
r — 2K/L, which destabilizes fluctuations about the 
ground state of the polymer. 

For an elastic filaments with torsional elastic con- 
stant T the total angle of rotation between two ends 
of a polymer is given by the sum of the writhing and 
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twisting contributions. The twisting contribution, 9 
is Gaussian distributed: 



Vjie) 



T 
27tl' 



=^T/2L 



(24) 



We can find the distribution of the total rotation an- 
gle cp from a convolution of the writhe and twist dis- 
tributions 

■p((p) = VK{yVr) Vj [Q] SilnWr + d - (p) dWr dd . 

(25) 
The twisting and writhing distributions are plotted 
in fig. (S). 




Figure 3: Distribution functions eq. (0), (EJ) of angle of 
rotation,x (on a log-linear scale) due to writhe and twist. 
The parabolic curve is the distribution of twist, the curve 
with the exponential wings is the distribution of rotation 
due to writhe. Plot is for T/L = K/L = 2 



5.3 Long Polymers 

The writhe properties of long semiflexible polymers 
have been extensively studied in the context of DNA 
super-coiling ||lq, Hq, n^. It is known from numer- 
ical simulations that for long semiflexible polymers 
(Wr^) ~ L. Recently, detailed analytic calculations 
| ]l8[ have shown that this is only true in models with 
a microscopic cut off. In the absence of a cut off 
the writhe fluctuations are logarithmically divergent. 
How do we understand these results from our picture 
of a stiff polymer as a random walk on the sphere? 

As shown by Fuller M the writhe of an arbitrary 
configuration of a polymer is found from the expres- 
sion 



Wr^-!- 

ZTT 



(to Xt) 
1 +t.to 



where to is an arbitrary constant vector. This ex- 
pression is valid when the path t(s) can be deformed 
from the vector to without the denominator vanish- 
ing. Let to point in the direction of the north pole 
then the singularity is at the south pole. To under- 
stand the nature of the singularity expand about the 
south pole so that 

t«(u,v,-1+uV2 + v2/2) (27) 

with both u and v small. The contribution to the 
writhe coming from a section of polymer near the 
south pole has the following form 



AWr 



1 



7t 



u dv — V du 1 



u^ 



7t 



dcf) 



(28) 



where cf) is the azimuthal angle of the path subtended 
at the south pole. The contribution of the polymer 
to the writhe near the pole is twice the winding num- 
ber about the pole. This winding number is clearly 
discontinuous as a function of the position of a path: 
A small circle centered at the pole has winding num- 
ber unity. If it is displaced more than its radius the 
winding number jumps to zero. 

We conclude that small fluctuations in the direc- 
tion of the polymer can give important (0(1 )) fluctu- 
ations in the total writhe if they occur near the south 
pole. The distribution of winding of a random walk 
on a sphere has been given analytically [EO, 19 . It is 



dt 



(26) 



known that the distribution is wide and has a Cauchy 
tail. This is due to the fact that in a close approach 
to the pole to within a distance e there are typically 
-yioge rotations about the pole. This is the origin of 
the Cauchy tail in the writhe fluctuations which de- 
cays even more slowly than the curve plotted in fig. 

(iO- 

Is this singularity at the south pole important for 
short polymers? When the polymer is very short 
transit in the region of the south pole is strongly sup- 
pressed since a bend of 180° costs an energy 7T^K/2L 
leading to exponentially small corrections to the law 
(pit). However as soon as L ~ K the dense nature of 
random walk in two dimensions means that with high 
probability there is a point on the polymer which is 
very close to the south pole, giving a strong Cauchy 
tail to the writhe distribution. 



6 Writhing dynamics 

When studying the dynamics of writhing of stiff poly- 
mer the representation of the writhe given by eq. (P6| 




Figure 4: Realization of a closed random walk on a 
sphere. The walk has a fractal nature with detail on all 
length scales. The area enclosed by the loop will be dom- 
inated by the largest loops in the above figure so that the 
area enclosed scales as the radius of gyration squared. 



is only useful for filaments which are much shorter 
than the persistence length. Motion of a long poly- 
mer is almost surely going to cause passage of the 
path over the south pole leading to a discontinuous 
contribution to the writhe. It is better to use the 
generalization ||] of eq. ( P6|) 



>Vr(t2)->Vr(ti) = -!- 



ti X t2 
1 +ti.t2 



(ti +i2] ds 



(29) 

Let the trajectory ti (s) to be the polymer at t = 
and t2(s) is the trajectory at some later time. This 
expression remains valid until a time such that a point 
on t(s) rotates 180° leading to a singularity in the 
integral. 

The motion of a semiflexible filament is given by 
the Langevin equation pi[ 



9r^ ^d'^r± 

at os^ 



+ f±{t,s) 



(30) 



This equation has a characteristic dispersion relation 
cu ~ q^ so that after a time t, wavelengths out to li — 
(Kt/r|)^/^ are equilibrated. Perturbations separated 
by more than the distance li have not had time to 



interact so that there are Ni — L/l] dynamically 
independent sections remaining on the polymer. 

What does this mean for the structure of the path 
t(s) on 52? From eq. ( po| ) the length li in real space 
translates to a distance Ri ~ yTi/K on 52. Structure 
on the scale of Ri has had time to move. We can thus 
imagine the structure at time t as being a series of 
N 1 beads of size -^/li/K strung out along the original 
path of the polymer on 52. The typical fluctuation 
in area from each bead is given by (AA) = ±R^. To 
find the total fluctuation in area in time t add these 
Ni independent fluctuations to find 



((>Vr(t)->Vr(0))^ 



Ni(AA)^ 
Lti/4 . 



Ui/k2(31) 
(32) 



This result has been found recently in detailed nu- 
merical simulations [p2| . The derivation in terms of 
random walks shows that the result is valid for poly- 
mers much longer than the persistence length: Indeed 
the result is valid for all times such that li (t) ^ K. 
Beyond this time there are two sources of breakdown 
in the theory: Firstly the dynamics are no-longer 
described by a simple bending theory and secondly 
we also expect that rotations by 1 80° become impor- 
tant leading to a breakdown of our expression for the 
writhe. 

Since the writhe fluctuations are the sum of a large 
number of independent contributions dynamic writhe 
fluctuations are distributed according to a Gaussian 
law. However there is always a exponentially small 
probability that the singularity in the integral occurs 
very rapidly so that the expressions found here is only 
be asymptotic, as was the case for the static distri- 
bution of writhe for short filaments eq. (E3) found 
above. 



7 Writhing of bent polymers 

The discussion until now has been restricted to poly- 
mers which are straight in their ground states. It 
is known, however, that certain special sequences of 
DNA lead to ground states which are strongly curved. 
This leads to enhanced writhing. Consider a poly- 
mer which bends an angle W over a distance L^, in 
its ground states. On the sphere the ground state is 
a geodesic of length W. We have seen above that 
a straight polymer is described by a diffusion like 
Fokker-Planck equation on the sphere. Addition of 
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drift due to the spontaneous curvature leads to the 
equation 



9P(t1 



W 



1 



= — V.VP + — V" 
9s Lb 2K ' 



.P(t1 



(33) 



with V a unit vector describing the direction of the 
bend. We take as our reference configuration ti the 
ground state and the state t2 the thermahzed poly- 
mer. As shown above the typical fluctuations of the 
filament on 52 scale as Rb ~ -y/Lb/K so that the 
area enclosed by the pair of paths ti and t2 scales 
as ±^Rb . Thus we find that 



(Wr^ 



K 



(34) 



This implies that a polymer, such as DNA with a 
sharp hairpin is expected to display enhanced tor- 
sional fluctuations at its ends due to this amplifica- 
tion of the writhe fluctuations by the spontaneous 
curvature. 

Similar modifications occur in the writhe dynam- 
ics: As before there are Ni = Lb/li dynamically in- 
dependent sections on a polymer, each section moves 
a distance a/Ii/K. It is stretched a distance ^li /Lb 
by the drift. Thus writhe fluctuations scale as 



{{Wr(t)-Wr(0)f 



- Ni((li¥/Lt 



)v/WK)f35) 
(36) 



showing that this effect is most important in short, 
sharp bends. 



8 Writhing of closed loops 

Recently simulations have been performed [Q to 
study the writhing dynamics of short closed DNA 
loops. Closure of the path in real space corresponds 
to adding a global constraint on the tangent 



t(a)^=0 

K 



(37) 



for all times. Here we consider the effect of this con- 
straint on the statics and dynamics of writhe fluctu- 
ations for a short polymer such that L ^ K. 

The zero temperature configuration of the polymer 
on 52 is a great circle which we take to be the equator. 
From eq. ( |37| ) the average vertical position on the 
sphere can not change: Any motion of the polymer 



in the north-south direction must be compensated by 
a corresponding south-north motion elsewhere in the 
polymer so that the average vertical position of the 
polymer does not move. However on the sphere this 
is a mass weighted average with the "mass density" 
1/k (Note: The enclosed area is calculated with a 
unit weighted measure). By coupling fluctuations in 
density along the equator with transverse fluctuations 
in the position of the equator we can find a non zero 
contribution to the writhing. 

The lowest order fluctuation in longitudinal den- 
sity which does not violate the constraint of eq. ( p7|) 
is proportional to cos2c() where cj) is the azimuthal 
angle. From eq. (pSJ ) its amplitude scales as a/L/K. 
Coupling this density fluctuation to the transverse 
motions, also of magnitude Rq- — ^L/K leads to a 
typical writhe fluctuation of ±L/K for a closed loop. 
Thus we find that 



(Wr^ 



lVk^. 



(38) 



How does the closure constraint modify the 
writhing dynamics? The geometry is very simi- 
lar to that discussed above for the bent open fila- 
ment. There are Ni = N/li dynamic sections of 
length ZttIi /L moving a distance -y/li/K. Again the 
first order contribution must vanish, since the "mass 
weighted" position of the equator can not move. A 
non-zero term can be generated by coupling to the 
longitudinal density fluctuations of amplitude -^/L/K: 



{[Writ] -Wr{0)f) ^l^/K^ 



h1/2 



(39) 



The absolute magnitude to the fluctuations is thus 
decreased compared with the open polymer but the 
dynamic exponent for the fluctuations remains 1/2. 
Recently [^ molecular dynamics simulations were 
performed on a detailed microscopic model of DNA 
loops. The numerical results seem to be consistent 
with a square root of time evolution of the writhe- 
writhe correlation function. 



9 Open Filaments 

9.1 Relation with geometric phases 

The geometry of fig. (|I|) is identical to that used to 
discuss the propagation of polarized light along bent 
fibre optics ||2^, ^, |2^ : Consider a circular cross sec- 
tion optical fiber bent into the shape shown in fig. 



9.1 Relation with geometric phases 



(mj. If plane polarized light is sent down the fiber 
with the plane of polarization defined by the vector 
n then the polarization will be transported in exactly 



the same manner as the normal vector to the bar 1 10 
After the triple bend of fig. (Q) the plane of rotation 
of the light has been rotated by 90°. The rotation 
of the plane of rotation in fibre optics has been ex- 
tensively studied as a simple example of the Berry or 
geometric phase. Berry's results on the existence of 
non-trivial phase factors is mathematically related to 
Fuller's theorem for writhe. While the geometry of 
the two situations seems very similar the geometric 
phase corresponding to this experimental situation is 
somewhat disguised so we shall now show how the 
rotation of the plane of polarization of light is equiv- 
alent to a phase shift for a photon. 

There are two equivalent mathematical descrip- 
tions available for polarized light: Firstly, planar po- 
larized waves, secondly circularly polarized photons. 
Using as a basis the planar representation with eigen- 
states 4)H — (1,0) and cfjv = (1,0), the helical eigen- 
states are (\)r — -7j(^,i) and 4)i, = 4=(1,— i). Both 
the pair {cjjv, cf^h} and the pair {(\>i, (\>r} form a com- 
plete orthonormal basis for describing arbitrary pure 
states of the light. Consider and experimental set up 
that rotates the vector (ph by an angle t|). The new 
state is given by (cosiJ;,sintj)). This state can also 
be expressed as a linear combination of the circularly 
polarized states. 



cosij; 
sintj) 



exp(-i\|;) 



1 



-exp(i\|;) 



1 



(40) 
The rotation of the plane of polarization is due to the 
modification of the relative phases of the two helical 
states. One state is advanced by ^\> the second is 
retarded ijj. When a beam of circularly polarized light 
is sent down a fiber there is a phase shift which is 
equal to the writhe of the fibre times the helicity of 
the photon. This is a Berry phase. 

A treatment of the geometry of writhe via a phase 
leads to an immediate generalization of Fuller's the- 
orem for trajectories which are not closed on the 
sphere [pm M: The rotation angle between n(0) and 
n(L) is for the moment undefined if t(0) and t(L) 
are not parallel. This is because we have only de- 
fined parallel locally on 52 via parallel transport. We 
shall now see how find a global definition of parallel 
by noting that the phase phase difference between 
two helical states is unambiguous even if the tra- 



jectory does not close on Sz'- The phase shift in a 
fiber with arbitrary boundary conditions corresponds 
to the following recipe for calculating the writhe of 
the filament (2|]: Follow the path t(0) to t(L) then 
close the path from t(L) to t(0) by a great circle (or 
geodesic). The phase shift is given by the area en- 
closed by the path t(0) --> t(L) augmented by the 
great circle t(L) — > t(0). 




v(L) 



Figure 5: Geometry of an open trajectory on Sz- 

The vectors t(0), t(L) and n(0) can be chosen to be co- 
planar. The vector v(L) is a reference vector at the point 
s = L which allows one to compare rotations with s = 0. 
The "most parallel" choice possible for v(L) is shown in 
the figure, and is coplanar with the other vectors in this 
figure. 

This recipe, fig. (^ has a interpretation on the 
sphere 52, for the writhing of a stiff polymer. Con- 
sider a path t(0) to t(L). The choice of n(0) is ar- 
bitrary for the calculation of the rotation, for conve- 
nience let us choose n(0) to that it points towards 
t(L). The three vectors t(0), t(L) and n(0) define a 
plane. Choose a reference vector, v(L), on the sphere 
at t(L) so that it is the "most parallel" possible to 
n(0). We choose v(L) so that it is coplanar with 
t(0) t(L) and n(0). This definition oi parallel ior two 
points on the sphere separated by a finite distance al- 
lows a definition of writhe for arbitrary free boundary 
conditions, leading to the geodesic recipe for calcu- 
lating the writhe. 
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9.2 Perturbative treatment of open 
filaments 

For short filaments there is an ahernative deriva- 
tion of the result that the path on the sphere must 
be closed by a geodesic. We shall use standard 
techniques from time dependent quantum mechanics, 
valid in the limit L/K <C 1 . We proceed by studying 
eq. (0) for the case of a general matrix, n(s), de- 
scribing both bend and twist, with Q^j — etjkUj'^. If 
the average direction of the polymer is aligned in the 
^z direction the angular velocities cu" and cu^ corre- 
spond to bending the filament, cu^ is the twist. We 
iteratively integrate this equation. The lowest order 
writhing contribution to rotations about the z-axis 
comes in second order. 



Q(Wr) 



dsds' 0(s 



(41) 
This is just the time ordered commutator of the ro- 
tation operators familiar from quantum perturbation 
theory. Integrating by parts to transform the 0- 
function into a 6-function gives 



a(Wr) 



'j^^ 



(tAt)ds-ht{0)At(L) (42) 



For a closed curve on 52 the first term is the area 
enclosed by the curve t(s) on a patch of a sphere in 
agreement with Fuller's theorem: It is the first term 
in the expansion of eq. (pq). The second boundary 
term, due to the integration, present when the curve 
is open corresponds to closing the path by a geodesic 
from t(L) to t(0). 

9.3 Dynamics of open filaments 

This definition of writhe has interesting consequences 
on the measurement of writhe fluctuations on freely 
fluctuating filaments. When the ends of the filament 
were constrained we showed that the writhing was 
due to a series of Ni small fiuctuations distributed 
along the chains. If t(L) can also fluctuate there is a 
new contribution due to the modification in the clos- 
ing of the path. In a time t t(L) moves yTi/K sweep- 

/L/K thus there is an extra 



ing out an area a/Ii/K x 

contribution to the writhing dynamics which scales 

(Wr^) ~ li L/K^ as before but coming from a single 

large end event rather than the sum of independent 

contributions. 



10 Conclusion 

We have shown that the writhing properties of a 
stiff polymer can be calculated from the well known 
properties of Gaussian random walks on a sphere. 
The predictions that (Wr^) - \J- and ((Wr(t) - 
>Vr(0))^) ~ t^/^L for short filaments should be ac- 
cessible to micromanipulation techniques on actin fil- 
aments, or perhaps DNA, either by studying the ro- 
tational dynamics of the end of a marked polymer 
or by direct observation of the shape fluctuations oc- 
curring in three dimensions via confocal microscopy. 
As pointed out in recent work Q disorder can some- 
times have strong effects on the dynamics of twist 
and torsional modes. It would be interesting to have 
a better understanding of the effect of disorder on the 
correlation functions discussed here. 

All real biopolymers have a helical structure with 
anisotropic bending elasticity. It is quite easy to treat 
this problem in the present description of diffusion 
on a sphere by introducing anisotropic diffusion coef- 
ficients on the sphere. 

I wish to thank P. Olmsted for discussions on the rela- 
tionship between geometric phases and Fuller's theorem, 
N. Rivier for discussions on Fermi- Walker (or parallel) 
transport, D. Beard for sending me his raw simulation on 
writhe fluctuations of DNA loops, A. Lesne for references 
on the distribution of winding angles. 
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